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Figure  1 .  Flow  cytometry  using  the  Life/Dead  assay  (Invitrogen)  for  single  species 
(a)  Detrick- 1  and  (b)  Detrick-2  cultures  after  24  h  of  incubation  in  1.5%  CSL.  The 
measured  fluorescence  intensity  from  SYT09  is  shown  in  the  y-axes  and  the  measured 
front  scattering  is  shown  in  the  x-axes.  The  left  graphs  represents  cells  stained  with 
SYT09  only  and  the  right  graphs  represents  cells  stained  SYT09  after  the  addition  of  PL 
The  red  colored  regions  correspond  to  the  live  cells,  the  blue  colored  regions  correspond 
to  the  dead  cells  or  the  cells  with  compromised  membrane  pumps,  and  the  green  regions 
correspond  to  the  cell  debris  or  other  particles . 6 

Figure  2.  Primer  set  quality  and  design  requirements,  (a)  Agarose  gel  of  the  products  from 
four  real-time  PCRs  using  Detrick- 1  and  Detrick-2  gDNA  as  template.  The  primer  sets 
corresponding  to  each  band  are  as  shown,  (b)  Traces  of  measured  fluorescence  intensity 
vs.  cycle  number  for  real-time  PCRs  using  Detrick- 1  and  Detrick-2  templates  and  one  of 
the  three  primer  sets  as  shown . 8 

Figure  3.  Sensitivity  of  the  PCRs  using  the  internal-control  primer  set  and  the  two  SNP- 
specific  primer  sets  and  gDNA  as  template  from  either  Detrick- 1  or  Detrick-2.  The 
highest  amount  of  template  per  PCR  is  10  ng,  which  corresponds  to  4,400,000  copies  per 
reaction.  The  template  copy  number  of  each  following  reaction  is  obtained  via  (1:9)  (v:v) 
serial  dilutions.  The  LightCycler480  system  was  used . 9 

Figure  4.  Specificity  of  primer  sets  pDl  and  pD2.  The  same  serial  template  dilutions  as  the 
ones  used  in  figure  3  are  used  here.  The  maximum  amount  of  mismatched  template 
undetectable  by  the  two  SNP-primer  sets  is  0.01  ng  or  less,  which  corresponds  to  4,400 
gDNA  copies  or  less . 10 

Figure  5.  Measuring  selectivity  of  SNP-specific  primer  sets  pDl  and  pD2.  Each  PCR 
reaction  contains  10  ng  total  of  gDNA;  however,  the  matching  template  from  either 
Detrick- 1  or  Detrick-2  is  different.  In  the  top  series,  the  following  ratios  of  Detrick- 
l:Detrick-2  (w:w)  templates  are  used:  (1:0),  (1:10),  (1:100),  (1:1,000),  (1:10,000), 
(1:100,000),  and  (1:1,000,000).  In  the  bottom  series,  the  same  ratios  are  used  for  the 
Detrick-2:Detrick-I  (v:v)  templates . 11 
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Figure  6.  Hypothetical  competition  experiments  for  the  emergence  of  (a)  Detrrick-1  and  (b) 
Detrick-2  in  mixed  cultures.  Left  graphs  compare  the  measured  values  from  mixed 

cultures  (red  data  points)  with  the  Cj  values  from  the  serial  dilution  experiments  in  figure 

5  (blue  data  points).  The  right  graphs  show  the  calculated  In  ( )  using  the  measured 

Cj  values  for  mixed  cultures  (equation  1)  vs.  hypothetical  generation  time.  The  ratios  of 

Detrick- 1  and  Detrick-2  cells  used  in  the  calculation  of  In  are  top  (Detrick- 

l:Detrick-2)  =  (0.002),  (0.01),  (0.02),  (0.1),  (0.5)  and  (1.0);  bottom  (Detrick-2:Detrick-l) 
=  (0.002),  (0.01),  (0.02),  (0.1),  (0.5)  and  (1.0).  The  slope  of  this  dependence  is  the  rate 
of  change  in  frequency  or  the  selection  coefficient  (for  non-interacting  unlimited 
cultures).  The  calculated  selection  coefficient  from  the  PCR  data  for  emerging  Detrick- 1 
is  equal  to  0.61 1  and  for  Detrick-2  it  is  equal  to  0.606.  The  theoretical  fitness  parameter 
(51)  is  0.694 . 
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1.  Objective 


Identification  and  detection  of  genomic  signatures  has  become  important  tool  in  studying  the 
development  of  drug  resistances  by  pathogens  or  the  emergence  of  high  morbidity  and/or  high 
mortality  pathogenic  strains.  The  ability  to  perform  competition  experiments  using  native  strains 
is  essential  for  understanding  the  conditions  for  strain  emergence  in  natural  environments  or  the 
interaction  of  pathogens  with  their  hosts.  This  report  develops  methodology  for  analyzing 
bacterial  competition  experiments  or  experiments  of  strain  emergence  using  native  strain- specific 
genomic  signatures.  The  goal  is  to  provide  a  sensitive  and  fast  method  for  measuring  relative 
fitness  parameters  of  non-engineered  bacterial  strains  with  sufficient  selectivity  to  allow  for 
environmental  applications.  The  proposed  methodology  has  potential  applications  in  forensic 
genomics,  where  the  contribution  of  a  given  mutation  to  strain  fitness  is  essential  for  inferring 
intent  in  markerless  engineered  genomes.  The  relative  polymerase  chain  reaction  (PCR) 
quantification  method  is  chosen  as  the  application  platform  for  the  proposed  methodology  and 
two  historically  important  military  strains  of  Bacillus  atrophaeus.  Detrick- 1  and  Detrick-2,  are 
used  as  an  example.  Detrick- 1  and  Detrick-2  are  congenic  strains  with  single  (C:T)  substitution 
corresponding  to  the  SpoOFHlOlR  mutation.  Our  accomplishments  include  the  development  of 
quantitative  assay  for  detection  of  the  (C:T)  single  nucleotide  polymorphism  in  spoOF,  which  is 
suitable  for  high-throughput  screening  of  mixed  bacterial  cultures.  The  assay  has  selectivity 
1:1000,  which  means  that  one  hypersporulating  bacteria  carrying  the  lOlR  signature  of  SpoOF  is 
detected  among  1000  bacteria  lacking  this  signature.  The  utility  of  the  assay  is  validated  in 
hypothetical  competition/emergence  experiments  for  either  Detrick- 1  or  Detrick-2.  The  assay 
has  been  applied  successfully  in  a  separate  work  for  studying  the  effect  of  SpoOHlOlR  mutation 
on  strain  fitness  (i). 


2.  Introduction 


B.  atrophaeus  var  globigii  (BG)  is  a  sporulating  soil  bacterium  used  as  a  common  bio-tracer  (2- 
6)  and  simulant  (7-15).  Recently  we  performed  whole  genome  molecular  typing  of  several  BG 
strains  including  the  strains  used  by  U.S.  military  as  simulants  for  biological  warfare  and 
bioterrorism  events  (16).  The  whole  genome  molecular  typing  of  the  military  strains  showed 
that  a  single  (C:T)  substitution  of  spoOF  corresponding  to  SpoOFHlOlR  mutation  supports  high 
spore  yield  (hypersporulating)  phenotype.  In  natural  environments,  sporulation  in  wild  type 
bacillus  species  is  commonly  induced  by  starvation  (17).  The  decision  to  commit  to  either 
sporulation  or  competence  is  a  complex  process  (18),  where  the  multicomponent  phosphorelay 
plays  a  key  role  (19)  and  where  SpoOF  is  a  major  regulator  for  phosphoryl  group  transfer  from 
several  histidine  kinases  to  SpoOB  (20)  and  then  to  the  transcription  regulator  SpoOA  (21),  which 
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regulates  the  expression  of  more  than  120  genes  (22).  The  SpoOF  protein  of  BG  is  identical  to 
the  same  protein  of  Bacillus  subtilis  except  for  two  amino  acids.  Similar  directed  mutation  at  the 
HlOl  residue  of  SpoOF  of  B.  subtilis  (SpoOFHlOlA)  confers  a  hypersporulating  phenotype  (23). 
It  is  believed  that  mutations  at  HlOl  alter  the  arrangement  of  the  amino-residues  in  the  p4-a4 
recognition  loop  and  affect  the  phosphorelay  performance  (24).  The  SpoOFHlOlA  mutation  in 
B.  subtilis  is  a  human-introduced  mutation.  In  contrast,  the  mutation  at  the  HIOIR  locus  of 
SpoOF  in  BG  occurred  spontaneously.  It  is  known  that  the  high  spore  yield  BG  strains  were 
intentionally  selected  in  experiments  conducted  after  the  Second  World  War  (25).  However,  it  is 
not  known  whether  these  strains  became  fixed  in  the  experimental  cultures  under  specific 
conditions  where  the  HlOlR  mutation  was  beneficiary.  The  inability  to  produce  a  viable  B. 
subtilis  strain  carrying  a  stable  HlOlR  mutation  suggests  that  this  mutation  might  not  be 
beneficiary,  in  which  case  the  fixation  under  natural  conditions  of  the  hypersporulating  strains  is 
unlikely  (16).  However,  the  inability  to  generate  a  viable  B.  subtilis  strain  carrying  the  HlOlR 
mutation  could  be  species  dependent. 

Bacteria  use  diverse  strategies  for  survival  in  changing  environments  including  genomic 
alterations  (26-28),  phenotypic  switching  (29),  or  mobilization  of  cooperative  responses  (30). 
Mutations  in  key  pathway  regulators  are  rare  events,  but  they  represent  an  important  mechanism 
for  increasing  phenotypic  diversity  (31).  Recent  directed  evolution  studies  combined  with  whole 
genome  sequencing  have  shown  that  relatively  small  number  of  mutations  in  key  regulators 
appearing  during  the  first  generations  under  new  selective  pressure  accounts  for  most  of  the 
gains  in  fitness,  and  that  subsequent  mutations  often  have  diminishing  effect  on  cumulative 
fitness  gains  (52-55).  The  acquisition  of  two  independent  mutations  in  SpoOF  of  BG,  where 
each  mutation  gives  rise  to  a  hypersporulating  variant  (16),  is  one  such  example  of  phenotype¬ 
defining  mutation  in  a  key  pathway  regulator.  One  approach  for  determining  whether  a 
phenotype-defining  mutation  is  beneficiary  in  specified  growth  conditions  is  to  compete  a  strain 
carrying  the  mutation  of  interest  with  its  congenic  ancestral  strain.  The  BG  pair  Detrick- 1  and 
Detrick-2  is  a  congenic  pair,  where  Detrick- 1  is  hyposporulating  strain  considered  ancestral  for 
the  hypersporulating  strain  Detrick-2.  Detrick-2  has  one  (C:T)  substitution  in  spoOF 
corresponding  to  the  SpoOFHlOlR  mutation  (16).  To  compete  the  two  strains,  their  frequencies 
in  the  competing  culture  need  to  be  measured.  Competition  experiments  using  antibiotic 
resistant  strains  have  been  performed  for  B.  subtilis  strains  obtained  from  directed  evolution 
experiments  (36,  37).  However,  the  use  of  antibiotic  resistance  could  introduce  unwanted 
complications  to  data  analysis  because  spores  are  naturally  resistant  to  antibiotics.  Therefore, 
here  single  nucleotide  polymorphisms  (SNP)-specific  quantitative  PCR  (qPCR)  assay  is 
developed  for  measuring  strain  frequencies  in  mixed  cultures,  which  accounts  for  the  cells  in  the 
culture  at  all  stages  of  their  life  cycle.  The  (T:C)  SNP  corresponding  to  the  HlOlR  mutation  is 
used  as  a  marker  for  distinguishing  between  the  two  strains.  SNP  detection  in  mixed  samples  is 
commonly  performed  using  either  end-point  detection  methods  (38-40)  or  real-time  PCR 
(qPCR)  (41-46).  Here  SNP-specific  qPCR  assay  based  on  the  relative  quantification  method 
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(47)  is  developed  and  used  for  calculation  of  the  rate  of  change  in  frequency  (48)  in  a 
hypothetical  experiments  of  emergence  of  either  Detrick- 1  or  Detrick-2. 


3.  Materials  and  Methods 


3.1  Bacterial  Cultures 

The  archival  strains  Detrick- 1  and  Detrick-2  of  BG  were  obtained  as  described  elsewhere  (16). 
Cells  were  cultured  in  corn  steep  liquor  (CSL)  (Sigma- Aldrich,  St.  Louis,  MO).  The  stock  CSL 
contained  approximately  50%  solids  and  50%  liquid  with  dissolved  water  soluble  components 
(49).  The  water  soluble  components  were  used  for  the  preparation  of  culture  media  as  described 
elsewhere  (25,  50).  Briefly,  the  stock  CSL  was  centrifuged  at  14,000xg  for  15  min.  The 
supernatant  was  collected  and  its  pH  was  adjusted  to  pH  7.0  with  sodium  hydroxide  (NaOH). 
Finally,  the  supernatant  was  cleared  by  performing  one  additional  centrifugation  at  14,000xg  for 
20  min.  The  cleared  supernatant  was  used  for  the  preparation  of  1.5  wt.%  CSL  growth  media. 
Alternatively,  the  cells  were  cultured  in  lysogeny  broth  (LB). 

To  avoid  possible  cross-contamination,  single  colonies  exhibiting  the  spoOF(lOlH)  or 
spoOF(lOlR)  signatures  were  chosen.  Three-day  or  older  colonies  were  used  for  starting  new 
cultures.  The  cells  from  a  single  colony  were  heat-shocked  and  grown  in  1.5%  CSL  as 
planktonic  cultures  to  optical  density  (OD)  between  1.0  and  2.0.  The  cultures  were  shaken  on  a 
rotary  shaker  at  250  rpm  and  33  °C  (MaxQ  4000,  Thermo  Scientific,  Asheville,  NC).  Then, 
0.5-ml  samples  were  taken  from  the  single-species  cultures  and  used  to  determine  the  percentage 
of  live  cells  using  Invitrogen’s  LiveDead  BacLight  assay  (Invitrogen,  Carlsbad,  CA)  and  flow 
cytometry  (BD  FACSCanto  II,  BD  Biosciences,  San  Jose,  CA).  Cultures  with  94%  or  higher 
percentage  of  live  cells  were  used  for  preparing  mixed  bacterial  cultures.  Strain  frequencies  in 
the  mixed  cultures  were  adjusted  using  OD  measurements. 

3.2  DNA  Isolation 

First,  6-ml  samples  from  either  single  species  cultures  or  mixed  cultures  were  collected  and 
washed  with  phosphate  buffered  saline  (PBS)  to  remove  cell  debris.  The  samples  were  washed 
five  additional  times  with  PBS  containing  10%  dimethyl  sulfoxide  (DMSO)  to  remove  possible 
polymerase  inhibitors.  Genomic  DNA  was  isolated  using  the  PowerSoil  kit  (MoBio 
Laboratories,  Carlsbad,  CA)  according  to  the  manufacturer  protocol  with  two  modifications. 

One,  10  pi  of  Proteinase  K  (Sigma)  10  mg/ml  was  added  to  solution  CL  The  cells  were 
incubated  in  solution  Cl  for  10  min  at  65  °C  prior  to  beating  in  a  bead  beater  for  40  s.  Second,  an 
additional  desalting  wash  with  ethanol  was  introduced  after  the  wash  with  solution  C5.  The 
quality  of  the  obtained  DNA  was  checked  using  NanoDrop  ND-1000  spectrophotometer 
(Thermo  Scientific). 
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3.3  Real-time  PCR 


SNP-specific  forward  primers  corresponding  to  C:T  substitution  of  spoOF  were  custom  designed 
using  the  sequence  information  from  the  whole  genome  typing  study  {16).  The  forward  primers 
were:  pSpo0F(101H)-Dl-f  matching  the  Detrick-1  DNA  strand 

(5’-TCAAAAGGCTTAGCGAAcT-3’)  and  pSpoOF(  10 lR)-D2-f  matching  the  Detrick-2  DNA 
strand  (5’-TCAAAAGGCTTAGCGAAcC-3’).  One  mismatch  in  -1  position  was  introduced 
(lower  case  nucleotide)  to  improve  primer  performance.  The  nucleotide  at  3’  position 
corresponded  to  the  (C:T)  SNP  of  the  spoOF  gene  {16).  One  common  reverse  primer 
pSpo0F(H101R)-186-rev  (5’-GCCCTGATCTTGTTCTCCTC-3’)  was  used  to  produce  186  bp 
amplicon.  Primers  specific  for  BG  and  producing  228-bp  amplicon  for  the  putative  RNA 
helicase  {16)  were  designed  and  used  as  internal  control  primers.  They  were  pBG21 1 15-f  (5’- 
AACCGTCGTAAACCAAGTCG-3’)  and  pBG21115-r  (5’-  TCCAAAGGTTCTGTCCCATC- 
3’).  Three  primer  sets  were  assembled  using  the  above  primers:  set  pDl  comprised  of 
(pSpoOF(lOlH)-Dl-f  and  pSpo0F(H101R)-186-rev  ),  set  pD2  comprised  of  (pSpo0F(101R)-D2-f  and 
pSpo0F(H101R)-l 86-rev  ),  and  set  p21115  comprised  of  (pBG211 15-f  and  pBG21115-r). 

PCR  reaction  volumes  were  10  pi,  including  500  pM  of  each  primer,  a  template  (template 
concentration  was  10  ng  if  not  specified),  and  5  pi  of  SsoFast-EvaGreen  supermix  (BioRad, 
Hercules  CA).  PCR  amplification  and  fluorescence  detection  were  carried  out  on  either 
LightCycler  (Roche)  or  CFX-Connect  (BioRad)  using  the  white  PCR  plates  for  LightCycler. 

The  thermal  cycling  for  LightCycler  included  a  5-min  initial  denaturation  at  94  °C, 
annealing/extension  for  40  s  at  54  °C,  and  fluorescence  acquisition  at  the  end  of  each  extension 
for  45  cycles.  The  thermal  cycling  for  CFX-Connect  included  a  5-min  initial  denaturation  at 
94  °C,  step  down  to  60  °C  for  3  s,  annealing/extension  for  40  s  at  54  °C,  and  fluorescence 
acquisition  at  the  end  of  each  extension  for  45  cycles.  Melting  curve  measurements  were 
performed  after  the  completion  of  cycling  of  each  real-time  PCR  reaction.  The  measured  Cj 
values  measured  using  either  the  LightCycler  or  the  CLX-Connect  systems  were  comparable. 


4.  Results 


4.1  SNP-specific  qPCR  Assay  for  Calculation  of  Fitness  Parameters 

In  competition  experiments,  two  bacterial  species  are  initially  mixed  at  a  desired  ratio  and  then 
incubated  in  specified  growth  conditions  for  several  generations  to  determine  whether  one  of  the 
species  has  fitness  advantage  in  the  chosen  conditions.  The  species  with  the  faster  growth  rate  is 
considered  to  have  a  fitness  advantage,  which  is  characterized  by  the  increase  of  the  species’ 
frequency.  The  logarithm  of  the  frequency  ratio  plotted  against  the  generation  time  is  used  to 
calculate  the  rate  of  change  in  frequency.  This  is  similar  to  the  calculation  of  selection 
coefficient  {51)  or  selection  rate  constant  (52)  of  non-interacting  cultures  during  the  exponential 
growth  phase.  Lor  density  dependent  cultures  or  cultures  grown  in  fluctuating  environments. 
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real  time  is  used  instead  of  generation  time.  The  logarithm  of  the  frequency  ratio  is  calculated 
from 

in  =  ln(l  +  (1) 

where  Px(.t)  is  frequency  of  strain  x  at  time  (t),  E  is  the  efficiency  of  the  PCR  reaction,  and 
^T,x{t)  C'T,y(t)  are  the  Cj  values  for  strains  x  and  y,  respectively.  By  convention  Detrick- 1  is 
considered  to  be  the  x  strain  and  Detrick-2  the  y  strain. 

4.2  Percent  Live  Cells  in  the  Initial  Cultures 

The  PowerSoil  kit  was  used  here  for  genomic  DNA  (gDNA)  isolation  from  both  cultures 
containing  both  vegetative  cells  and  spores.  The  Invitrogen’s  Live/Dead  assay  was  used  to 
measure  the  percentage  of  live  cells  in  the  single- species  starting  cultures  of  Detrick- 1  and 
Detrick-2  (figure  1).  The  left  graphs  in  figure  1  show  the  measured  fluorescence  from  SYT09 
versus  front  scattering  for  cells  strained  with  SYT09  only.  The  right  graphs  in  figure  1  show  the 
measured  SYT09  fluorescence  after  the  addition  of  propidium  iodide  (PI). 
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Figure  1.  Flow  cytometry  using  the  Life/Dead  assay  (Invitrogen)  for  single  species  (a)  Detrick- 1  and 

(b)  Detrick-2  cultures  after  24  h  of  incubation  in  1.5%  CSL.  The  measured  fluorescence  intensity 
from  SYT09  is  shown  in  the  y-axes  and  the  measured  front  scattering  is  shown  in  the  x-axes.  The 
left  graphs  represents  cells  stained  with  SYT09  only  and  the  right  graphs  represents  cells  stained 
S  YT09  after  the  addition  of  PI.  The  red  colored  regions  correspond  to  the  live  cells,  the  blue 
colored  regions  correspond  to  the  dead  cells  or  the  cells  with  compromised  membrane  pumps,  and 
the  green  regions  correspond  to  the  cell  debris  or  other  particles. 

SYT09  is  a  minor-groove  binding  fluorophore,  which  has  high  membrane  permeability  and 
stains  both  live  and  dead  cells.  PI  is  another  minor  groove  binder,  which  quenches  the 
fluorescence  from  SYT09.  PI  is  pumped  out  of  the  cells  when  the  membrane  pumps  are 
functional,  which  prevents  the  intracellular  PI  concentration  to  reach  levels  sufficient  for 
quenching  SYT09  fluorescence.  This  effect  is  used  to  distinguish  between  live  cells  or  cells 
with  functional  membrane  pumps  and  dead  cells  or  cells  with  compromised  membrane  pumps. 
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Three  distinct  regions  are  identified  in  the  graphs  shown  in  figure  1 :  (a)  red-colored  region 
where  SYT09  fluorescence  is  not  affected  by  the  addition  of  PI  identifying  the  live  cell 
population;  (b)  blue-colored  region  where  the  addition  of  PI  quenches  the  SYT09  fluorescence 
identifying  the  dead  cell  population  or  the  cells  with  compromised  membrane  pumps;  and 
(c)  green-colored  region  where  the  SYT09  fluorescence  is  at  low  level  before  the  addition  of  PI 
identifying  cell  debris  or  other  particles.  The  data  in  figure  1  show  that  96.3%  of  Detrick- 1  cells 
and  92.8%  percent  of  Detrick-2  cells  are  live  cells,  respectively.  The  dead  cell  population 
represents  2.4%  of  the  Detrick- 1  population  and  4.3%  of  the  Detrick-2  population,  respectively. 
These  data  show  that  more  Detrick-2  than  Detrick- 1  cells  are  either  dead  cells  or  have  non¬ 
functional  membrane  pumps.  This  percentage  of  live  cells  was  taken  under  consideration  during 
the  construction  of  the  mixed  cell  cultures  for  the  hypothetical  emerging  experiments  discussed 
below. 

4.3  Requirements  for  Primer  Sets  Design 

The  efficacy  of  the  relative  quantification  PCR  method  is  sensitive  to  multiple  factors  including 
polymerase  used,  primer  quality,  reaction  conditions,  template-dependent  bias,  amplicon  size, 
etc.  The  utility  of  the  method  is  dependent  on  the  ability  to  minimize  the  differential  effect  of 
these  factors  when  SNP-specific  primer  sets  are  used.  Figure  2a  shows  the  agarose  gel  of  the 
products  generated  in  four  PCRs.  Two  of  the  reactions  use  gDNA  from  Detrick- 1  as  template 
and  two  use  gDNA  from  Detrick-2  as  template,  respectively.  Two  primer  sets  are  used  in  each 
reaction;  one  is  a  SNP-specific  primer  set  and  the  other  is  the  internal  control  primer  set.  Two 
bands  are  seen  in  the  lanes  where  the  SNP-specific  primer  sets  are  matching  the  used  templates. 
One  of  the  bands  corresponds  to  the  SNP-specific  primer  set  and  the  other  corresponds  to  the 
internal-control  primer  set.  Only  one  band  is  seen  in  the  lanes  where  the  SNP-specific  primer 
sets  are  mismatching  the  used  template.  This  band  corresponds  to  the  internal-control  prime  set. 
No  other  bands  appear  on  the  gel,  showing  that  the  designed  primer  sets  are  suitable  for  use  in 
real-time  PCR  assays. 

The  SNP-specific  primer  sets  are  designed  to  fulfill  the  requirement  for  minimum  bias  toward 
the  used  templates.  This  is  illustrated  in  figure  2b,  where  the  Cj  for  the  matching  templates  of 
the  two  SNP-specific  primer  sets  are  almost  equal.  The  traces  from  the  real-time  PCRs  using  the 
internal-control  primers  are  also  shown  as  reference.  The  internal-control  primer  set  is  designed 
choosing  a  region  in  the  BG  genome  with  specified  GC  content  and  size,  such  that  the  measured 
Cj  values  when  the  internal-control  primer  set  is  used  to  be  similar  to  the  Cj  values  when  the 
SNP-specific  primer  sets  are  used  with  matching  template.  This  provided  for  the  opportunity  to 
use  directly  the  measured  Cj  values  for  the  reactions  using  SNP-specific  primer  sets  in  equation 
1  without  adjustments  for  template  bias. 
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A 


tDl  tDl  tD2  tD2 

pDl  pD2  pDl  pD2 


Figure  2.  Primer  set  quality  and  design  requirements,  (a)  Agarose  gel  of  the  products  from  four  real-time 
PCRs  using  Detrick- 1  and  Detrick-2  gDNA  as  template.  The  primer  sets  corresponding  to  each 
band  are  as  shown,  (b)  Traces  of  measured  fluorescence  intensity  vs.  cycle  number  for  real-time 
PCRs  using  Detrick- 1  and  Detrick-2  templates  and  one  of  the  three  primer  sets  as  shown. 

4.4  Primer  Sets  Sensitivity,  Specificity,  and  Selectivity 

The  performance  of  SNP-specific  primer  sets  is  characterized  by  sensitivity,  specificity,  and 
selectivity.  Sensitivity  is  defined  as  the  minimum  number  of  template  molecules  detected  by  the 
assay.  The  theoretical  limit  for  PCR  sensitivity  is  one  template.  This  limit  is  easily  approached 
in  the  reactions  using  the  internal-control  primer  set  p21 1 15  when  the  reactions  using  the  SNP- 
specific  primer  sets  are  characterized  by  lower  sensitivity. 

Figure  3  shows  the  results  from  experiments  used  to  determine  assay  sensitivity.  The  reactions 
with  the  highest  Cj  values  shown  in  figure  3  use  10  ng  of  gDNA  per  reaction,  which  corresponds 
to  4,400,000  copies  of  gDNA.  The  template  copy-number  in  each  flowing  reaction  is  obtained 
via  (1:9)  (v:v)  serial  dilutions  using  PCR-clean  water.  The  lowest  amount  of  template  per 
reaction  is  obtained  after  six  serial  dilutions,  which  corresponds  to  four  templates  per  reaction. 
Figure  2a  and  b  shows  that  when  the  p21 1 15  primer  set  is  used  the  dependence  of  the  measured 
Cj  values  on  the  log  of  template  concentration  is  linear  over  the  full  range  of  serial  dilutions 
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used.  In  contrast,  when  the  SNP-specific  primer  sets  are  used  the  above  linear  dependence  is 
observed  only  over  the  range  of  template  copies  from  4,000,000  to  40. 


Template  -  tDl;  Primer  set  -  p21115 


Template  -  tD2;  Primer  set  -  p21115 


Template  -  tDl;  Primer  set  -  pDl 


Template  -  tD2;  Primer  set  -  pD2 


Figure  3.  Sensitivity  of  the  PCRs  using  the  internal-control  primer  set  and  the  two  SNP-specific  primer  sets 
and  gDNA  as  template  from  either  Detrick- 1  or  Detrick-2.  The  highest  amount  of  template  per 
PCR  is  10  ng,  which  corresponds  to  4,400,000  copies  per  reaction.  The  template  copy  number  of 
each  following  reaction  is  obtained  via  (1;9)  (v:v)  serial  dilutions.  The  LightCycler480  system 
was  used. 
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One  additional  requirement  is  introduced  for  all  reactions,  which  is  to  have  efficiency  close  to 
100%.  The  efficiencies  of  the  reactions  using  the  internal-control  primer  set  were  93.7%  for  the 
Detrick- 1  template  and  95.2%  for  the  Detrick-2  template.  The  reaction  efficiencies  increased 
and  approached  100%  when  the  SNP-specific  primer  sets  are  used.  Therefore,  the  additional 
requirement  is  fulfilled  by  all  primer  sets. 

Specificity  is  measured  for  reactions  using  SNP-specific  primers  and  is  defined  as  the  maximum 
amount  of  mismatched  template  for  which  the  measured  Cj  values  are  equal  to  the  Cj  values  for 
the  control.  Figure  4  shows  the  results  of  experiments  performed  in  the  determination  of  primer 
set  specificity.  The  specificity  of  the  pDl  primer  set  is  on  the  order  of  4,400  copies  and  the 
specificity  of  the  pD2  primer  set  is  on  the  order  of  2,200  copies  or  less. 

Template -tDl  Template -tD2 

Primer  set  -  pD2  Primer  set  -  pDl 


Figure  4.  Specificity  of  primer  sets  pDl  and  pD2.  The  same  serial  template  dilutions  as  the  ones  used  in 
figure  3  are  used  here.  The  maximum  amount  of  mismatched  template  undetectable  by  the  two 
SNP-primer  sets  is  0.01  ng  or  less,  which  corresponds  to  4,400  gDNA  copies  or  less. 

Selectivity  characterizes  the  ability  of  SNP-specific  primer  sets  to  respond  to  the  simultaneous 
presence  of  matching  and  mismatching  templates.  Selectivity  is  coupled  to  specificity;  however, 
in  practice,  it  is  defined  via  the  ratio  of  the  frequency  of  the  matching  and  mismatching 
templates.  Alternatively,  selectivity  is  expressed  as  the  percentage  of  detectable  matching 
templates.  Here  selectivity  is  measured  as  the  ratio  of  matching  and  mismatching  templates. 

The  measured  selectivity  for  the  two  SNP-specific  primer  sets  shown  in  figure  5  is  at  least 
(1:1000).  This  means  that  one  copy  of  the  matching  template  is  detected  in  the  presence  of  1000 
copies  of  the  mismatching  template.  From  figure  5,  the  selectivity  of  primer  set  pDl  is  1:1000 
and  the  selectivity  of  primer  set  pD2  is  1:4000.  The  efficiency  of  the  reactions  within  the  range 
of  selective  template  detection  is  1 10%  for  primer  set  pDl  and  102.8%  for  primer  set  pD2. 
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Template  -  tDl  in  tD2;  Primer  set  -  pDl 


Efficiency 

2.100 

Slope 

-3.103 

Intercept 

20.10 

Enoi 

0.123 

Template  -  tD2  in  tDl;  Primer  set  -  pD2 


Efficiency 

2.024 

Slope 

-3.266 

Intercept 

20.31 

Enor 

0.119 

Figure  5.  Measuring  selectivity  of  SNP-specific  primer  sets  pDl  and  pD2.  Each  PCR  reaction  contains  10  ng 
total  of  gDNA;  however,  the  matching  template  from  either  Detrick- 1  or  Detrick-2  is  different.  In  the 
top  series,  the  following  ratios  of  Detrick- 1: Detrick-2  (w:w)  templates  are  used:  (1:0),  (1:10),  (1:100), 
(1:1,000),  (1:10,000),  (1:100,000),  and  (1:1,000,000).  In  the  bottom  series,  the  same  ratios  are  used 
for  the  Detrick-2:Detrick-l  (v:v)  templates. 

4.5  Analysis  of  Competition/Emerging  Experiments 

The  application  of  the  relative  quantification  PCR  method  for  analysis  of  competition/emerging 
experiments  is  illustrated  using  a  hypothetical  emerging  experiment,  where  Detrick- 1  emerges  in 
a  mixed  culture  with  Detrick-2  or  vice  versa.  The  strain  with  the  smaller  initial  frequency  is 
considered  to  be  the  “emerging”  strain.  When  Detrick-2  is  the  “emerging”  strain  the  experiment 
is  labeled  as  “Detrick- 1  elimination,”  because  by  convention  Detrick-2  is  the  y  strain  for  the 
purposes  of  using  equation  1 .  One  unit  of  hypothetical  generation  time  is  defined  as  the  time 
necessary  for  the  doubling  of  the  emerging  strain.  At  generation  time  to,  the  frequency  of  the 
emerging  strain  is  set  equal  to  0.001.  Using  these  conditions,  the  calculated  theoretical  selection 
coefficient  (51)  is  equal  to  0.693. 

The  mixed  cultures  shown  in  figure  6  are  constructed  to  mimic  the  conditions  for  strain 
emergence  in  mixed  cultures.  OD  measurements  were  used  during  the  construction  of  the 
genotypically  mixed  cultures  used  in  the  experiments.  gDNA  was  isolated  from  the  mixed 
cultures  and  the  total  DNA  concentration  was  adjusted  to  10  ng  per  reaction  using  the  internal- 
control  primer  set  p21 1 15.  The  Cj  values  for  the  reactions  using  the  matching  SNP-specific 
primer  sets  were  measured  and  placed  together  with  the  Cj  values  for  the  serial  dilution 
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experiments  in  figure  5.  The  species  mixing  ratios  (see  legend  of  figure  6)  were  used  to 
determine  the  position  of  the  measured  Cj  values  on  the  x-axes.  The  combined  data  sets  are 
shown  as  left  graphs  in  figure  6.  The  graphs  show  that  the  Cj  values  measured  for  the  mixed 
cultures  are  in  very  good  agreement  with  their  counterparts  measured  for  the  template  dilution 

series.  Therefore,  the  efficiencies  from  figure  5  are  used  for  the  calculation  of  In  ( )  in  the 

right  graphs  in  figure  6.  The  selection  coefficients  calculated  from  the  slopes  of  the  dependence 
of  In  ( )  on  the  hypothetical  generation  time  shown  in  figure  6  were  equal  to  0.61 1  for  the 

hypothetical  emergence  of  Detrick- 1  and  0.606  for  the  hypothetical  emergence  of  Detrick-2. 
These  values  are  close  to  the  theoretical  selection  coefficient  equal  to  0.694  calculated  above. 
This  demonstrates  the  utility  of  the  proposed  relative  PCR  quantification  approach  for  analysis  of 
competition  experiments. 


A  tDl -f  tD2  =  20  ng  per  10  ul  reaction 


hypothetical  emergence  of  Detrick-1 


log(tDl/tD2) 


g  tDl  +  tD2  =  20  ng  per  10  ul  reaction 


log{tD2/tDl) 


hypothetical  elimination  of  Detrick-1 


y 

=  -0.6064X  +  7.2834 

=  0.991 

theor.  fitness  param.  = 

-0.694 

4  6  8  10 

hypothetical  generation  time 
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Figure  6.  Hypothetical  competition  experiments  for  the  emergence  of  (a)  Detrrick-1  and  (b)  Detrick-2  in  mixed 
cultures.  Left  graphs  compare  the  measured  values  from  mixed  cultures  (red  data  points)  with  the 
values  from  the  serial  dilution  experiments  in  figure  5  (blue  data  points).  The  right  graphs  show  the 

calculated  In  { using  the  measured  Cj-  values  for  mixed  cultures  (equation  1)  vs.  hypothetical 
\Py(£)/ 

generation  time.  The  ratios  of  Detrick-1  and  Detrick-2  cells  used  in  the  calculation  of  In  ( )  are  top 

VPy(£)/ 

(Detrick-1  ;Deti-ick-2)  =  (0.002),  (0.01),  (0.02),  (0.1),  (0.5)  and  (1.0);  bottom  (Deti-ick-2:Detrick-l)  = 
(0.002),  (0.01),  (0.02),  (0.1),  (0.5)  and  (1.0).  The  slope  of  this  dependence  is  the  rate  of  change  in 
frequency  or  the  selection  coefficient  (for  non-interacting  unlimited  cultures).  The  calculated  selection 
coefficient  from  the  PCR  data  for  emerging  Detrick- 1  is  equal  to  0.61 1  and  for  Detrick-2  it  is  equal  to 
0.606.  The  theoretical  fitness  parameter  (57)  is  0.694. 
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5.  Conclusions 


A  novel  approach  for  analysis  of  competition  experiments  based  on  the  use  of  SNP-specific 
primer  sets  and  the  relative  PCR  quantification  method  is  developed.  The  utility  of  the  approach 
is  validated  using  the  BG  species  Detrick-2  and  Detrick- 1.  These  species  are  suitable  model 
system  for  single-locus  haploid  selection  experiments  because  they  are  isogenic  except  for  a 
single  (C:T)  SNP  corresponding  to  the  HIOIR  mutation  of  SpoOF  {16).  The  developed  assay 
detects  the  (C:T)  SNP  of  SpoOF  with  (1:1000)  selectivity,  which  is  sufficient  for  analysis  of 
competition  experiments,  or  the  detection  of  the  Army  simulant  strain  Detrick-2  or  its  progenitor 
species  in  environmental  samples. 
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BG 

Ct 

CSL 

Detrick- 1 
Detrick- 2 
DMSO 
E 

gDNA 

LB 

OD 

PAt) 

PBS 

PCR 

pDl,  pD2  and  p21115 
PI 

qPCR 

SNP 

SpoOA 

SpoOB 

SpoOF 

t 

tDl  and  tD2 


Bacillus  atrophaeus  var  globigii 

cycle-threshold 

corn  steep  liquor 

low  spore  yield  variant  of  Bacillus  atrophaeus  var  globigii 

high  spore  yield  variant  of  Bacillus  atrophaeus  var  globigii 

dimethyl  sulfoxide 

efficiency  of  real-time  PCR  reaction 

genomic  DNA 

lysogeny  broth 

optical  density 

frequency  of  a  species  in  mixed  cultures 

phosphate  buffered  saline 

polymerase  chain  reaction 

primer  sets 

propidium  iodide 

quantitative  PCR 

single  nucleotide  polymorphism 

transcription  factor  (response)  regulator 

phosphotransferase 

two  component  response  regulator 

time  measured  in  generations 

genomic  DNA  templates 
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